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We theoretically study the development of spontaneous magnetization in diluted magnetic semi- 
conductors as arising from a percolation of bound magnetic polarons. Within the framework of 
a generalized percolation theory we derive analytic expressions for the Curie temperature and the 
magnetization, obtaining excellent quantitative agreement with Monte Carlo simulation results and 
good qualitative agreement with experimental results. 
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Diluted magnetic semiconductors (e. g. Ini_2:Mna;As, 
Gai_a;Mna;As), which are materials with some fraction 
of the nonmagnetic lattice atoms replaced by magnetic 
atoms {i. e. a fraction x of Ga or In atoms being ran- 
domly replaced by Mn atoms), have attracted a great 
deal of attention 0, || |[ iTS, |, 0, |, |, fol- 
lowing the recent discovery ]11| of ferromagnetism in 
Gai-xMn^jAs, with a Curie temperature Tc « lOOK in 
the X = 0.03 — 0.07 range. The subject is of considerable 
fundamental and technological interest. Technologically, 
a semiconductor, which is also a ferromagnet, raises the 
exciting potential of spintronic applications where 
logic and memory operations could in principle be seam- 
lessly integrated on a single device. From a fundamen- 
tal perspective, understanding ferromagnetism in a novel 
material (which is also a semiconductor) is an important 
challenge, particularly so in this case in light of the ex- 
tensive existing earlier work involving doped II- VI semi- 
conductor materials (as against the systems of current 
interest, which are magnetically doped III-V systems) 
which were never convincingly found to be ferromagnetic 
except perhaps with extremely low values of Tc. It is 
therefore not surprising that there has been a great deal 
of theoretical activity |, | |, |, | |, | |, |lO| trying 
to understand the ferromagnetic mechanism in GaMnAs. 
Although no theoretical consensus on the precise ferro- 
magnetic mechanism has yet been reached in the liter- 
ature it is now established that the interaction between 
the magnetic atoms, which leads to the ferromagnetic 
phase at low enough temperatures, is induced by charge 
carriers (holes in the case of GaMnAs) in the semicon- 
ductor host. Theoretical investigation of these systems 
is hampered by the fact that both disorder and inter- 
actions are strong and must be taken into account non- 
perturbatively. This problem has been approached in a 
number of ways, including theoretical approximations as- 
suming charge carriers being almost free ^, |^, Q and 
numerical studies in the opposite limit of strongly local- 
ized charge carriers ^, ^ 0]. However, a comprehensive 
understanding of the physics of diluted magnetic semi- 
conductors has not been achieved yet. 

Theoretical approaches 0, ||, Q treating the charge 
carriers as free carriers in the valence band of the semi- 



conductor have employed the simple Weiss mean field 
theory incorporating the band structure details {e.g. 
spin-orbit coupling, etc.), without taking into account 
effects of disorder. While this approach is claimed to 
provide quantitatively accurate values of Tc, it fails qual- 
itatively in accounting for experimentally observed trans- 
port properties ||l3| of GaMnAs where the resistivity is 
very high and seems to obey a Mott variable range hop- 
ping behavior at low temperatures reminiscent of an in- 
sulating system. In addition, the resistivity always de- 
creases with increasing temperature above Tc, which is 
again typical of a localized insulating system. The free- 
carrier mean-field approach also fails to explain the re- 
cently reported dependence of Tc on annealing which 
indicates a crucial role for disorder. A recent dynamical 
mean field theory calculation concluded that the fer- 
romagnetic GaMnAs may actually belong to a situation 
where the carriers are just on the edge of being strongly 
localized. The approach ||, ^, 0| treating the carriers 
as localized carriers in a semiconductor impurity band 
(thus being opposite to the the free-carrier mean-field 
approach), which has so far been explored only numer- 
ically, also leads to reasonable agreement with experi- 
ments and indicates a very strong dependence of Tc on 
disorder. In this paper we provide an analytical the- 
ory which takes into account both disorder and strong 
magnetic interaction starting with localized carriers and 
using the physically appealing magnetic polaron perco- 
lation picture. Our theory has a starting point similar 
to that in Refs. |[ Q except that ours is a completely 
analytical physical theory in contrast to the numerical 
approach used in the works |6| 0]. Where applica- 
ble our analytical results agree well with the numerical 
results of Refs. [||, 0. 

In this paper we consider a system in which transi- 
tion to the insulating state due to localization of the 
charge carriers occurs at temperatures higher than the 
Curie temperature Tc. The carriers will be called "holes" 
throughout the text of this paper since in GaMnAs the 
carries are holes although the theory to be developed 
in this paper applies equally to the situation where the 
carriers are electrons. Exchange interaction of localized 
holes with magnetic impurities leads to the formation of 



bound magnetic polarons [T^. Since the concentra- 
tion of magnetic impurities is much larger than the hole 
concentration most likely due to compensation by 
As antisite defects, a bound magnetic polaron consists of 
one localized hole, and a large number of magnetic impu- 
rities around the hole localization center. Even though 
the direct exchange interaction of the localized holes is 
antiferromagnetic, the interaction between bound mag- 
netic polarons may be ferromagnetic at large enough 
concentrations of magnetic impurities. To understand 
the physics of this phenomenon, one may consider two 
neighboring polarons (Fig. |l|). The localized holes of 
these polarons both act on the impurities surrounding 
them thus producing an effective magnetic field for these 
impurities. The energy minimum is reached by this sys- 
tem when the impurity spins are parallel to this effec- 
tive field, and the magnitude of the field is maximum. 
The maximum of this effective magnetic field is achieved 
when the spins of the localized holes are parallel. There- 
fore at low temperatures the system should eventually 
reach the state where the spins of all holes point in the 
same direction, and all impurity spins point in the same 
or in the opposite direction, depending on the sign of the 
impurity-hole exchange interaction. 

If the hole localization radius is much less than the 
characteristic distance between the localized holes, the 
disorder in the hole positions must have dramatic ef- 
fect on the whole picture of the ferromagnetic transition. 
This notion has been confirmed by Berciu and Bhatt §, 
who have shown, by means of numerical simulations, that 
both the Curie temperature and the shape of the mag- 
netization curve M (T) are indeed strongly affected by 
disorder. It has been known that the percolation theory 
p8| provides many adequate tools to deal with ferromag- 
netism in disordered systems with strong localization of 
carriers . In this paper we present a quantitative de- 
scription of the spontaneous magnetization in magnetic 
semiconductors within the framework of the percolation 
theory. 

In our model, the charge carriers are localized. The 
hole wave function is assumed to fall off exponentially 
away from localization centers, with decay length as- We 
consider the low carrier density regime in which the mean 
distance between the localized holes is much larger than 
the hole localization radius, a%nh ^ 1. We note that for 
GaMnAs, as ~ 10 A, and therefore our theory applies in 
the regime rih <C 10^^ cm~^, with experimental rih values 
currently being around 10^^ cm^^. The localization cen- 
ters are distributed randomly in the sample. Magnetic 
impurities are distributed within the sample with con- 
centration Hi 3> rih randomly as well. The Hamiltonian 
of the system has the form 




kj 



(1) 



where indices k and j label magnetic impurities and holes 



FIG. 1: Interaction of two bound magnetic polarons (after 
Ref. [Q) The polarons are shown with gray circles, small and 
large arrows show impurity and hole spins respectively. 



respectively, Sk/sj are the impurity/hole spin operators. 
Matrix elements Jkj of the impurity-hole exchange in- 
teraction decay exponentially with the distance between 
the interacting impurity and hole, Jkj = Joexp(— 2|rfe — 
Tjl/as)- The direct (antiferromagnetic) exchange in- 
teraction between the magnetic impurities is neglected, 
since their relative concentration x in the lattice of the 
host semiconductor is much less than unity although this 
may be important for larger values of x, leading to the 
suppression of Tc. We will not discuss possible mech- 
anisms of hole localization, since the properties of our 
model hold for any of them as long as the decay of the 
localized hole's wave function is exponential. 

At some temperature T, magnetic impurities that are 
at distances r < Rp{T) = (05/2) ln(sS'| Jo|/T) from hole 
localization centers have their spins strongly correlated 
with the spins of the corresponding holes (here s and S 
are the absolute values of the hole or impurity spin re- 
spectively). The spins which do not have a localized hole 
within a circle of radius Rp(T) around them are essen- 
tially free. The quantity Rp{T) is the effective radius of a 
magnetic polaron; it grows as the temperature is lowered. 
Clearly at low enough temperatures neighboring mag- 
netic polarons overlap and interact with each other via 
interaction with impurities between them |p^ . This in- 
teraction produces alignment of the polaron spins. When 
the cluster of correlated polarons having the size of the 
sample (the infinite cluster) appears, the ferromagnetic 
transition occurs. 

Even though the picture presented above is qualita- 
tively correct, one needs to find the maximum character- 
istic temperature T2p{r) at which the spins of two po- 
larons at a given distance r from each other are still 
strongly correlated. The Hamiltonian of a two-polaron 
subsystem is given by Eq. ([T]), where hole index j takes 
only two values ji and j2 corresponding to the two po- 
larons under consideration. Our present goal is to find 
the characteristic spin-correlation temperature T2p{r) at 
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given r, as, Jo, and n^. Since we are interested in the 
system behavior only at and below the percolation tran- 
sition, we may limit our consideration to polaron pairs 



Hamiltonian to the following form: 



with 



> 



-1/3 



^ as- A rough preliminary estimate, 



which can be obtained from straightforward dimensional 
analysis, reads 



T2p{r) ^ A{r,aB,ni)sS\Jo\exp{-r/aB) 



(2) 



where dependence of the dimensionless prefactor func- 
tion A{r,aB,ni) on its arguments is weaker than expo- 
nential. Since at T ~ SiSj JqI exp(— r/a^) the two po- 
larons must be correlated, A > 1. At temperatures of 
the order of the right-hand side of Eq. (||), each polaron 
already has a large number of impurity spins near its 
center strongly polarized along the direction of the spin 
of its hole. Due to this coupling, the spin of the po- 
laron's hole becomes "massive," and may be considered 
a classical field as far as its interaction with more re- 
mote impurities is concerned. This effective field pro- 
duces polarization of the remote impurities; the charac- 
teristic value of the projection 5'^"'^^ of the spin of fcth 
impurity onto the spin direction of the ji-th polaron is 
roughly 5min{l, sJkji/T}. The characteristic energy of 
interaction of this impurity with the other polaron is of 
the order of sSJkj2, therefore the contribution of the 
fcth impurity to the interaction of two polarons is of the 
order of 5*^5 J^j^ min{l, sJ^j^/T} [here we assume that 
sSJkji,T > sSJkj2 since T2p(r) > sS\Jo\exp{-r/aB)]- 
This estimate allows us to determine which impurities 
are important for the interaction of two given polarons. 
Namely, these impurities are the ones in the lens-shaped 
region between the two localization centers, see Fig. ||. 
The diameter of this region is of the order of y'ras; more 
remote impurities are interacting too weakly with both 
polarons to be of any importance. Finding the thickness 
of the lens-shaped region is less trivial. As we consider 
impurities away from the middle point and closer, say, 
to ji, decrease of Jkj2 is fully compensated by the in- 
crease of Jkji , until we enter the region where the spins 
of the impurities are saturated due to their proximity to 
ji. Thus the size of the region of interest in the direc- 
tion along — r^^ depends on the temperature. Since 
the characteristic correlation temperature T2p(r) is the 
quantity to find, we will proceed by iterations, starting 
with the value given by the right-hand side of Eq. (|^) 
with A{r, aB,ni) = 1 (so the thickness of the lens-shaped 
interaction region equals as). As we will see, just one 
iteration yields the answer with good precision. 

We will take into account only the impurities in the 
interaction region; the total number of these impurities 
equals N ~ a^rrii in the first iteration. The coupling of 
each of these impurities to either hole is approximately 
J(r/2) = Jo exp(— r/as). It allows us to reduce the 



N 



H = sJ{r/2) COS -Y^Si'^ 



(3) 



k=l 



Here the angle between the spins of the two polarons 
is denoted 9, and the direction of the z axis is chosen 
along the direction of the vector Sj^ + Sj^ . We have ne- 
glected spatial variation of the polaron exchange field in 
the interaction region, since the relative magnitude of 
this variation does not exceed unity. 

The partition function of the system described by 
Hamiltonian (|^) can easily be found, and after straight- 
forward algebra we arrive at the expression for the av- 
erage cosine of the angle 9 between the spins of the two 
interacting polarons: 



N[sSJ{r/2)/TY 



1 - 



N[sSJir/2)/T]^ 



NsS\J{r/2)\ <r 
NsS\J{r/2)\ >r 



(4) 



Therefore, the spins of the two polarons at distance r are 
correlated at temperatures below y/NsS\J{r/2)\. Using 
the estimate N « a^rrii, we finally arrive at 



T2p{r) ^ aBy/rnisS\Jo\ exp (-r/aB) 



(5) 



The maximum possible distance VcorriT) between two 
polarons with correlated spins at a given temperature T 
is therefore 



r{T) ~ as 



In 



T 



1 



■In 



jn,; In 



sS\Jo\ 
T 



(6) 



Now we can use result (||) to get a better estimate for 
the thickness of the lens-shaped interaction region in or- 
der to perform the next iteration. It is determined by 
the condition Jkj^ ^ < T, so the thickness is of the order 
of as In (as y^rnj). The resulting correction to T2p{r) is 
just a factor of In(aB^rni). Taking this correction into 
account is clearly beyond the accuracy limits of the ap- 
proximation used above, so we neglect it, and use Eq. ^ 
as our final result. 

Now we consider a system of randomly placed magnetic 
polarons. As the temperature is being lowered, the spins 
of neighboring polarons become aligned, and clusters of 
polarons with the same spin appear. At any given tem- 
perature T, the polarons separated by a distance smaller 
than TcorriT) are joined into magnetic clusters. The 
lower the temperature, the more such "links" between po- 
larons are established, and the larger the average cluster 
size. Finally, at low enough temperatures, a cluster hav- 
ing the dimensions of the sample, the so called "infinite 
cluster" appears, and the magnetization of the sample ac- 
quires some finite value. The problem of finding the tran- 
sition temperature is identical to the problem of finding 
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the critical percolation radius in the problem of randomly 
placed overlapping spheres iQ. The latter problem has 
been solved numerically, and it has been demonstrated 
that the infinite cluster forms in the system when the link 
length reaches the value Tperc ~ 0.86/ ^/uh Substi- 
tuting this expression into Eq. (^), we get the expression 
for the ferromagnetic transition temperature 



T'c ~ asy^n^^^'^sS'l Jol exp 



0.86 

aB^/nh 



(7) 



The limit of applicability of Eq. (j^) is determined by the 
condition a'^n^ <C 1. A similar exponent was obtained 
in Ref. ||2^ for ferromagnetic transition in PdFe alloys. 

It is instructive to point out that our Eq. is consis- 
tent with the mean-field result derived in the literature 
0, H, ^, Q for the opposite limit of almost free holes, 



mf 



Jo, and 



where J^f is related to Jq of Eq. (jl|) by Jmf 
£,{T) ~ nil/ max{eF,T} is the magnetic susceptibility of 
holes, with sp being the Fermi energy. In the case of low 
electron density, T ^ ep 



Tc ^ VTVnh\Jmf\ , 

which matches the result (0) of the bound-polaron pic- 
ture at the limit of applicability of the latter, a^rih 1. 

Since Uh <C Ui, one polaron includes many magnetic 
impurities, and the total magnetization of the sample is 
that of impurities: 

M{T) = SniV[rcorriT,ni,nh,aB)\/TV^ , (8) 

where rcorr is defined by Eq. (H), the and universal func- 
tion V{y) is the infinite cluster's volume in the model 
of overlapping spheres; it depends only on the product 
y of the spheres' diameter and the cubic root of their 
concentration. 

Using Eqs. and (0), we cast Eq. ^ in the form 



M(T) 
M(0) 



V 0.86 -f (a|nh)Mn 



(9) 



with the Curie temperature Tc given by Eq. (^. One can 
see that the shape of the magnetization curve is deter- 
mined by only one dimensionless parameter a^n/j, while 
the expression for Tc is more complicated and de- 
pends on all parameters of the model. Fig. || shows the 
temperature dependence of the magnetization at two val- 
ues of a^rih] the curve is more concave at smaller values 
of this parameter, which is precisely the experimental 
observation. Our magnetization results agree with the 
numerical results of Ref. and are consistent with the 
experimental magnetization data in GaMnAs, particu- 
larly for systems with lower values of Tc where our po- 
laron percolation picture applies better due to stronger 
carrier localization associated with lower values of a^rih- 




FIG. 2: The solid lines show the relative magnetization of 
the magnetic impurities [Eq. (^] for a%nh = 5 x 10"'^ (cmve 
1) and 2 x 10~^ (curve 2). The dashed lines show the relative 
magnetization of localized holes, whose contribution to the 
total sample magnetization is small. 



To conclude we have developed an analytic polaron 
percolation theory for DMS ferromagnetism in the limit 
of low carrier density or equivalently strong carrier local- 
ization a'grih ^ 1. Interestingly our polaron percolation 
theory reproduces the free carrier Weiss mean field the- 
ory in the limit of a^rih ~ 1. Our analytic results are 
in good agreement with existing numerical results in the 
strongly localized limit |§, ^, . The experimental DMS 
currently have a^n/i ~ 10~^ — 10~^ which makes our the- 
ory marginally applicable to the experimental systems, 
and we get reasonable agreement with experimental re- 
suhs for Tc and for M(T). 
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